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Abstract 

The formulation of the low-momentum effective interaction in the model space Lee-Suzuki and 
the renormalization group methods is implemented in the three-dimensional approach. In this 
approach the low-momentum effective interaction Vi ow k has been formulated as a function of the 
magnitude of momentum vectors and the angle between them. As an application the spin-isospin 
independent Malfliet-Tjon potential has been used into the model space Lee-Suzuki method and it 
has been shown that the low-momentum effective interaction Vi ow k reproduces the same two-body 
observables obtained by the bare potential Vnn- 
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I. INTRODUCTION 



One of the most fundamental and essential problems in the nuclear theory is the derivation 
of the effective interaction between the nucleons of nuclei. The conventional nuclear force 
models are given by the same one-pion exchange interaction in the long distance, but differ 
substantially in their treatment of the short distance. The fact that the different short 
distance structures reproduce the same observables for two-body problem indicates that 
low energy observables are insensitive to the details of the short distance dynamics. This 
insensitivity is a result of the separation of the high and low energy scales in the nuclear 
force, and implies that we can derive a low-momentum effective interaction. 

During the past years, several methods have been developed to derive the energy inde- 
pendent low-momentum effective interaction such as the renormalization group (RG) and 
the model space techniques. These approaches are mainly based on a partial wave (PW) 
decomposition and the details have been given in references [1-11]. Bogner et al. have 
developed a low-momentum effective interaction which is quite successful in describing the 
two-nucleon system at low energy. This effective interaction is independent of the potential 
models as the cutoff is lowered to A = 2.1 fm~ x [l, 2|. 

Recently the three-dimensional (3D) approach, which greatly simplifies the numerical 
calculations of few-body systems without performing the PW decomposition, is developed 
for few-body bound and scattering problems [12-27]. The motivation for developing this 
approach is introducing a direct solution of the integral equations avoiding the very involved 
angular momentum algebra occurring for the permutations, transformations and especially 
for the three-body forces. Conceptually the 3D formalism consider all partial wave channels 
automatically. Based on this non PW method the bound and scattering amplitudes are 
formulated as function of momentum vector variables, specially their magnitudes and the 
angle between them. 

In this article our aim is to generate the low-momentum effective interaction directly in 
a 3D formalism, where as a simplification the spin and isospin degrees of freedom have been 
neglected in the first step. So the low-momentum effective interaction has been formulated 
for spinless identical particles as function of vector Jacobi momenta. Considering the spin 
and isospin degrees of freedom is a major additional task, which will increase more degrees of 
freedom into the states [25( . The formulation of the low- momentum effective interaction in a 
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spin-isospin dependent 3D approach, based on helicity representation, is currently underway 



and it will be reported elsewhere 



28|]. 



This article is organized as follows. In the section [IT] the model space Lee-Suzuki method 
has been used to derive the energy independent model space effective interaction in the 3D 
representation. In addition a renormalization group decimation method has been performed 
in the 3D representation and a flow equation for the low-momentum effective interaction 
Vi ow k has been obtained as a function of momentum vectors in the appendix El The section 
III 1 1 describes the details of the discussion and numerical calculations of the low-momentum 
effective interaction Vi ow k in the model space Lee-Suzuki method by using the spin-isospin 
independent toy Malfliet-Tjon potential. Finally a summary is given in the section [IV] and 
an outlook is provided. 



II. MODEL SPACE LEE-SUZUKI METHOD FOR V iowk IN THE 3D MOMENTUM 
REPRESENTATION 

Several model space techniques such as the Bloch- Horowitz [l|, Q); the Unitary Transfor- 
mation [3j, |5| and the Lee-Suzuki [7-11] based on the restricted space, have been developed 
for the construction of the low-momentum effective interaction Vi ow k- 

In this work the Lee-Suzuki formalism has been applied to the free space nucleon-nucleon 
problem in the 3D momentum representation and the low-momentum effective interaction 
Vi ow k has been obtained as function of the momentum vectors. A momentum-space Hamil- 
tonian for the full-space two-body problem has been considered as follows: 

H(k, k') = H (k)5(k - k') + V NN (k, k'), (1) 

where Ho denotes the kinetic energy and VjvAr(k, k') is the bare two-body interaction. In 
the model space methods the projection operators onto the physically important low-energy 
model space, the P space, and the high-energy complement, the Q space, have been intro- 
duced as: 



P = j rf 3 A;|k)(k|,|k| < A, 

Q = J d 3 k |k)(k|, |k| > A, (2) 



3 



where A is a momentum cutoff which divides the Hilbert space into the low and high mo- 
mentum states. These projections satisfy the relations: 



P + Q = 1, 
PQ = QP = 0, 
P 2 = P, 

Q 2 = Q, (3) 
and they act on the full-space two-body problem states as: 

p\^ N ) = i*r>, 

Q\i>n = (4) 

where ) and l^j^) denote the states of the full and model spaces respectively. The u> 
is an operator which transforms the states of the P space to the states of the Q space. The 
key aspect of the Lee-Suzuki method is the determination of the u> operator defined by the 



following equation 



29|: 



QM")=Q«>PM")- (5) 



The Schrodinger equation for the full-space two-body problem by considering the equation 
(JSJ) can be written as: 

H\^ N }=H(P + QM N } = H(P + QcoPM n }=E\^ n ). (6) 

By acting the P operator on the left side of the last equation the full-space two-body problem 
has been reduced to the model space two-body problem of the following form: 

PH(P + QujP)\^ n ) = (PH P + PVP + PVQujP)\^ n ) = EP\^ N ), (7) 

where we assumed that the P and Q operators commute with the Hq. Therefore the non- 
hermitian low-momentum effective potential in the model space that reproduces the model 
space components of the wave function from the full-space wave function can be written as 
[7-11]: 

V lowk = PV(P + QuoP). (8) 
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By using the integral form of the projection operators P and Q, the low-momentum effective 
interaction Vi ow k can be written in the 3D representation as: 

/OO P 
dqq 2 J dqVi vw (k , ,q)a;(q,k), (9) 

where q is the momentum vector in the complement model space Q. For an application we 
can choose the suitable coordinate system which the vector k is along the z axis and the 
vector k' is in x — z plane. So the equation ([9]) can be rewritten as: 



p OO /*1 /*27T 

Viowkik', k, x) = V NN (k', k,x) + / dqq 2 dx'j dip V NN (k',q,y) u(q,k,x'), 

Ja J -i Jo 



(10) 



where: 



X = 


k- k. 


x' = 


k- q, 


y = 


k • q 


k = 


|k|, 


k' = 


in 


q = 


1 q|- 



X 12 COS(f, 



(11) 

To calculate the low-momentum effective potential Vi ow k we need to determine to. To this 
aim by applying (q| on the left side of equation (JSJ) and using the integral form of projection 
operators P and Q, this equation can be rewritten as follows: 

*f(q) = £ dpp 2 J rfpc(q, p) ¥£"(p), (12) 

where p is the momentum vector in the model space. We use the completeness relation in 
the model space as: 

£dkk 2 Jdk V» N (p , )V™(p) = 5(p'-p). (13) 
The equation (fT2l after implementing equation (1T3"]) can be written as: 

w(q, p) = £ dkk 2 J dk ^(q)^(p), (14) 
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where ^^(p) and ^/^(q) are the wave function components of the P and Q spaces of the 
full-space respectively which in the form of the half-on-shell (HOS) two-body T matrix are 
given by: 



<"(P) = <5(P"k)+ r(p ' k ' P) 



k 2 — p 2 + is 
T( P ,k,^ 



= 8(p - k) + p ' ' - ztt - p 2 ) T(p, k, ^), (16) 
where p denotes a principle value. The HOS two-body T matrix can be obtained from the 



Lippmann-Schwinger equation in the 3D representation, which is given as 

7^*.) = W > + f^Y^mfl. (17) 

In order to solve the equation (1141) we have chosen suitable coordinate system which the 
vector q is along the z axis and the vector p is in x — z plane. Therefore the equation (fl4"|) 
can be rewritten as: 



u(q,p,x)= [ dkk 2 [ dx' r d^ N (q,x')^ N (p,y), (18) 
Jo J -i Jo 



where: 



X = 


q - P, 


1 

X = 


q - k, 


y = 


k P 


k = 


|k|, 


p = 


IpI> 


q = 


1 q|- 



x' 2 COS If, 



(19) 



We should mention that the Lee-Suzuki method reproduces the HOS two-body T matrix 
the same as the RG method, therefore the solution o : 
be equivalent to the solution of the RG equation 
renormalization group decimation method in the 3D representation and we obtain a flow 
equation for the low-momentum effective interaction V[ ow k for future applications. 



this approach has been proven to 

n n 

|4j, |30|]. In appendix (jAj) we derive a 
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FIG. 1: The total two-body cross section from the low- momentum effective potential (circles) and 
the bare potential (solid line) as a function of kinetic energy in the lab frame. 

III. DISCUSSION AND NUMERICAL RESULTS 



In our calculations we employ the spin-isospin independent Malfliet-Tjon potential. This 
force is a superposition of a short-range repulsive and long-range attractive Yukawa inter- 
actions. It is given as [31] : 

1 , V R V a 



Wk',k) 



(20) 



2tt 2 v (k'-k) 2 + /4 (k'-k) 2 + /4 7 ' 
With the strengths and the masses of the meson exchange as follow: Vr = 7.291, Va = 
3.177, /xr = 613.7MeV, /jla = 305.9MeV . This potential supports one bound state at E — 
— 2.23MeV. With this interaction we first implement the LU factorization into equation 
(fT3|) to calculate Of% N (j/, x', ip') as an inverse of the wave function ^^ N (p, x, <p) in the model 
space. Then by solving the equation (fT8|) we calculate u>(q,p,x) and finally we input the uj 
operator into equation ffTUl) to obtain the low- momentum effective interaction Vi ow k{p' ,p, x). 

The dependence on the continuous momentum and the angle variables is replaced in the 
numerical treatment by a dependence on the certain discrete values. To this aim we use the 
Gaussian quadrature grid points to discrete the momentum and the angle variables. In our 
calculations we choose forty grid points for the momentum variables in the P space, i.e. the 
interval [0, A], and thirty grid points for the momentum variables in the Q space, i.e. the 
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FIG. 2: The two-body 5-wave phase shift from the low-momentum effective potential (circles) and 
the bare potential (solid line) as a function of kinetic energy in the lab frame. 



interval [A, oo). Also twenty and fourteen grid points for the spherical and the polar angle 
variables have been used respectively. The integration interval for the P and Q spaces are 
covered by two different linear and tangential mappings the Gauss-Legendre points x from 
the interval (-1,+!) via: 



P 



A, 



q = btan (-(1 + x) J + A , b = 5 fm' 1 . 



(21) 



to the intervals [0, A] and [A, oo) respectively. As we mentioned in the introduction we have 
used the value of 2.1 fm^ 1 for the cutoff A in our calculations. The solution of the integral 
equation (TPTj) requires a one-dimensional interpolation on T. We use the cubic hermitian 
splines of ref. 32| for its accuracy and high computational speed. It can be useful to mention 



33(, for solving a system of 



that in the numerical calculations we use the Lapack library 
linear equations in the calculation of the two-body T matrix. 

We have shown in figure [T] that the calculated total two-body cross section for the 
Viowkip'jP, x ) an d the Vjvjv(p')P) x ) match perfectly well. Also we have calculated the S- 
wave phase shift for Vi ow k(p',p, x) and it has been compared with the obtained result with 
the bare potential Vnn(p',P, x) in figure [2j For calculation of the phase shift we have used 



S 



P [fm _1 ] o 1 

(a) 



P [fm 1 ] o 1 

(b) 



FIG. 3: (a) The comparison of the low-momentum effective potential Vi ow k(p-P, x) (solid lines) 
with the bare potential Vnn(p,P, x) (dashed lines) and (b) differences between them. 




FIG. 4: The comparison of uj(q,p)(q 2 — p 2 ) operator for the S'-wave calculated in the PW approach 
(dotes) and in the 3D approach (solid lines). 



the relation between the PW and 3D representations of the on-shell T matrix as follow: 

Ti(p,p,p 2 ) = 2tt J dx T(p,p,x,p 2 ) Pi(x). (22) 
As we expect the results are in good agreement with the high accuracy. In figure ([3^) we 



FIG. 5: The comparison of the low-momentum effective potential V lowk (p' ,p) for the S'-wave cal- 
culated in the PW approach (dotes) and in the 3D approach (solid lines). 

have compared the calculated low-momentum effective potential Vi ow k{p,p, x) with the bare 
potential Vnn(p,P, %), also In figure ([3b) the difference between them is shown as a function 
of the momentum variable p and the angle variable x. This comparison demonstrates that the 
difference between the two potentials is approximately constant except for the momentum 
values close to the cutoff, where the differences vary for different angle values. Conceptually 
when we integrate on the angles we can produce the same results as PW approach. 

As a test of our calculations we can compare the obtained results for the low-momentum 
effective potential and the u operator in both 3D and PW approaches. In the first step we 
directly calculate the u operator and the low-momentum effective potential in PW approach. 
In the second step we obtain a certain PW projection of the u operator and the low- 
momentum effective potential from their corresponding 3D representation by the following 
relations: 



The obtained results for the S'-wave projection of the u operator and the low-momentum 
effective potential in both of 3D and PW approaches are given in figures (J1J and fl5j) respec- 




(23) 



fO 



tively. The agreement between the two approaches is quiet satisfactory. 

IV. SUMMARY AND OUTLOOK 

In this article the 3D formulation of the Lee-Suzuki and the RG methods have been 
presented. The low-momentum effective interaction Vi ow k has been derived as a function 
of the magnitude of momenta and the angle between them without using the partial wave 
decomposition. The calculated two-body observables from the low-momentum effective in- 
teraction and the bare interaction have been shown. In addition a comparison of calculated 
Vi ow k from the PW and the 3D representation have been demonstrated as a test for our 
calculations. 

The advantage of our formulation in the 3D representation in comparison with the PW 
representation is that we have calculated the low-momentum effective interaction by consid- 
ering all partial waves automatically. 

For the future investigations the low-momentum effective interaction can be formulated 
by considering the spin and isospin degrees of freedom in a realistic 3D approach. This 
formulation based on the momentum helicity basis has been done and the calculation of 
the realistic low-momentum effective interaction with Bonn-B, AV18 and Chiral potentials 
is currently underway. Considering the obtained 3D low-momentum effective interaction in 
the few-body bound and scattering calculations is another major task to be done. 

Acknowledgments 

We would like to thank S. K. Bogner for fruitful discussion during EFB19 conference 
about this work. This work was supported by the research council of the University of 
Tehran. 

APPENDIX A: RG METHOD FOR V iowk IN THE 3D MOMENTUM REPRE- 
SENTATION 

For the two-body problem, the RG method derives the low-momentum effective potential 
Viowk by integrating out the model dependent high momentum component of the different 
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models of bare potentials Vat at by demanding that the low-momentum observables calculated 
from Vnn must be reproduced by V\ owk with the same accuracy We follow basically the 
same procedure which is used by Bogner et al. in derivation of the RG flow equation 4\. 
However instead of using the PW representation, we have implemented the treatment of the 
RG in the 3D representation. We start with a RG treatment of the scattering problem in 
the 3D representation. The HOS two-body T matrix in 3D representation by applying an 
arbitrary potentia. model^ to the scattering problem is given by the Lippmann-Schwinger 
equation (h = m N = 1) [12] [34]: 

T(k,k,fc ) = V NN (k,k) + pj dpp J dp k 2 - p 2 ' ^ ' 

where p denotes a principle value integration. Corresponding to equation (|7j) we can define 
a restricted version of the equation ( 1A1I) by imposing a cutoff A and replacing the bare 
potential Vnn by a low-momentum effective potential Vi ow k as: 

'' Viow k(k', p) T lowk (p, k, k 2 ) 
k 2 — p 2 

(A2) 



T lowk (k.', k, k 2 ) = Vi owk (k', k) + p dpp 2 J d£>- 



We demand that the calculated low- momentum HOS two-body T matrices from equations 
(fATl) and $55} are identical, i.e. T(k',k,£; 2 ) = T lowk {k', k, k 2 ) for |k'|, |k| < A. This 
condition ensures that the V\ owk gives the same low-momentum two-body observables as 
obtained by the bare potential Vnn- Imposing the cutoff independence of the low- momentum 
HOS two-body T matrix, i.e. dT(k', k, k 2 )/dA = in equation flA2j) . ensures that the 
low-momentum observables will be independent of the scale A and will generate energy- 
independent potential Vi owk , therefore we obtain: 

A 2 [ ^ dV lowk (k',p) f ^ V lowk (k',Ap)T(Ap,k,k 2 ) 
dpp / dp 77 Xk(p) = / dp , (A3) 



./ dA J r l-(k/A) 2 

where we have used the standing wave scattering states of the effective theory \xk), which 
can be written as: 

lx k ) = |k) + p^ A ^ P 2 | dp T( ^ k 'p |p>. (A4) 



[1] This quantity is real in this appendix 
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By using the completeness relation of the scattering states in the model space and considering 
T(Ap,k, k 2 ) = (Ap|VJ ow fc|xk), we have obtained: 

^(k'lV^felp') = A 2 J dp (k'|V^|Ap> 

= A 2 J dp(k'\V lowk \Ap) £ dp"p" 2 J dp"{Ap\V lowk \p") 



x Idkk 2 (d&"^fW 
Jo J A 2 - k 2 

= A 2 J dp(k'|V^ fe |Ap) 

x J Q A dp"p" 2 f d p" (Ap| 14™ fe I p") C(p", p', A 2 ) , (A5) 

where (xic| is the bi-orthogonal complement in completeness relation in the model space 
which satisfies (Xk'IXk) = S(k' — k) and G denotes the interacting Green's function in the 
model space. With writing the G in terms of the two-body T matrix we obtain the RG 
equation in the 3D momentum representation as: 

d M (V ^ [ ^ Vio W k(k',Ap)T(Ap,k,A 2 ) 

^Wk,k) = J dp f - wx? . (A6) 

For numerical calculations we can choose the suitable coordinate system where the vector k' 
is along the z axis and the vector k is in x — z plane. By this consideration we can rewrite 
the above equation as: 
d 



dA 

where: 



v (Vi, \ f 1 a i f 2 \ Vio W k(k',A,x>)T(A,k,y,A 2 ) 
V lowk {k ,k,x) = J dx j dip 1 _ (k/A) 2 ' ^ ^ 



x — k * k ; 



X 



k'p, 

y = k • p = xx' + \f\ — x 2 \/l — x' 2 cosip, 
k = |k|, 

k' = \k'\. (A8) 



We introduce: 



f(k',k,x',x) — / dipT(k',k,x'x + Vl-x' 2 Vl-x 2 cosip). (A9) 
Jo 
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And finally the RG equation in the 3D representation can be obtained: 

- ,*,*) = ^ rfx f - WXj , . (A10) 

The low- momentum effective interaction Vi ow k(k',k,x) can be calculated by numerically 
integrating the RG equation on the Gaussian momentum and angle grid points. Similar to 
the PW representation of this equation, VNN(k',k,x) can be used as a large-cutoff initial 
condition to calculate Vi ow k{k', k, x) numerically. 
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